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We investigate theoretically, the character of electronic eigenstates and transmission properties 
of a one dimensional array of stubs with Cantor geometry. Within the framework of real space 
re-normalization group (RSRG) and transfer matrix methods we analyze the resonant transmission 
and extended wave-functions in a Cantor array of stubs, which lack translational order. Apart from 
resonant states with high transmittance we unravel a whole family of wave-functions supported 
by such an array clamped between two-infinite ordered leads, which have an extended character 
in the RSRG scheme, but, for such states the transmission coefficient across the lead-sample-lead 
structure decays following a power-law as the system grows in size. This feature is explained from 
renormalization group ideas and may lead to the possibility of trapping of electronic, optical or 
acoustic waves in such hierarchical geometries. 

PACS numbers: 73.21. La,73.63.Kv,85.35.Be 



I. INTRODUCTION 

Electron transmission in quantum wires with various 
geometric structures has been an interestingarea of re- 
search in mesoscopic physics in recent times [li Q, y, |j, |M 
0, 0, IE • Apart from presenting an interesting problem 
in itself, this area of research has profound implications 
in applied physics and materials science, particularly, in 
designing nano-devices with exotic electronic properties. 

Nano-technology, as its present advanced stage, makes it 
possible to fabricate a variety of quantum structures such 
as quantum dots, wires and rings. The surface of solids 
can be designed with appreciable control to serve as ideal 
templates for the formation of low-dimensional nano- 
structures with desired geometry |0 . Nano-devices such 
as one-dimensional wires can be experimentally prepared 
by controlling self-organized defects at the surface [lll |. 
Apart from the path-breaking experiments, from the 
point of view of a theoretician, an interesting observa- 
tion is that, often the basic principle underlying the op- 
eration of these fabricated nano-systems can be modeled 
in a quite simple manner, using quantum mechanics and 
elementary solid state physics. One such example is the 
study of electron transport in model quantum wires with 
dangling side branches (stubs) 2., Ai A, 3» Jji> M of 
finite length. Such 'dangling' stubs can be grafted along 
a nano-wire with any specified ordering. The presence 
of stubs in a quantum wire has been shown to lead to 
transistor action by Sols et al 0. The length of the 
stubs can be controlled by applying a gate voltage, and 
a change in these lengths result in a non-trivial change 
in the spectrum of the system 2, 3]. The spectrum and 
band structure of an array of stubs were studied by Deo 
and Jayannavar as well. They pointed out that, un- 
like a system of potential scatterers, a system of geomet- 



rical scatterers (such as an array of stubs) exhibit com- 
plete band formation in the conductance with only a few 
such stubs. Several studies on the quantum waveguide 
transport using model stubbed nano- wires then followed. 
Shi and Gu j5( examined the electronic conductance of 
a quantum waveguide with side-branches taking into ac- 
count the impedance factors for geometric and potential 
scatterers. Vasseur et al 0] studied simple tight bind- 
ing models of one dimensional comb structures of simple 
metals. In a similar spirit Pouthier and Girardet It'I have 
investigated the conductance of a monatomic nano-wire 
containing a random distribution of side-grafted nano- 
clusters of atoms. Cervero and Rodriguez modeled 
a quantum wire by a random array of delta potentials 
and extended their work to investigate the effect of posi- 
tionally correlated potentials on the spectral behaviour. 
Orellana et al 9] stuck to the tight hinging model to 
investigate the electronic transport through a quantum 
wire with a side quantum dot array. Quite recently, 
Cao et al 0] have studied the electronic transport in 
a comb structure with a mesoscopic ring threaded by a 
magnetic flux. While the Schrodinger equation approach 
0, 0, IE Hi 0j Q is a natural choice for the calculation 
of the conductance, discrete difference equation approach 
lESSl also finds its justification by the recent advance- 
ment in atom manipulation technology in which one can 
place individual atoms on the substrate by using a scan- 
ning tunnel microscope tip as tweezers @. 



Though the transmission characteristics of a stubbed 
quantum wire are quite sensitive to the geometry of the 
structure, little effort has been given to study the spectral 
characteristics when the arrangement of the stubs devi- 
ates from being periodic. One interesting area, where 
in-depth studies in this regard are really lacking is, when 
the stubs are arranged in a hierarchical pattern. In this 
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communication we undertake a detailed study of the elec- 
tronic states and transport in a model one-dimensional 
quantum wire with stubs of finite length grafted on it 
followin g a triadic Cantor sequence (a one-dimensional 
fractal) |l3| . Apart from providing direct experimental 
evidence for localized and self-similar fracton mode dis- 
placements interesting theoretical results indicate 
the presence of multiple fracton and phonon regimes in 
a triadic Cantor sequence 0|. Also, the Cantor distri- 
bution of dielectric layers has recently found application 
in the design of guided-wave optical or microwave filters 
|15|. Comparatively speaking, the electronic properties 
of a Cantor quantum waveguide network is unattended so 
far, though recent tight binding results on the electronic 
properties of a Cantor lattice with point scatterers in- 
dicate exotic features being markedly different in 
several aspects from the conventional quasi-periodic lat- 
tices of the Fibonacci class . In the case of a stubbed 
waveguide network the possibility of tuning the stub- 
lengths experimentally may provide an opportunity for 
practically 'seeing' some of the peculiar features already 
encountered in the theoretical studies of such Cantor-like 
substitutional sequences [TgIIT^ . 

With this motivation we undertake a detailed analysis of 
electronic transmission and the nature of wave-function 
in a Cantor stub waveguide (CSW) network using the 
real space renormalization group (RSRG) and the trans- 
fer matrix (TM) method. We discretize the Schrodinger 
equation to map the original sequence of stubs onto 
an equivalent one-dimensional chain consisting of point- 
like scatterers placed at the nodes where a stub meets 
the backbone. This facilitates the application of RSRG 
method. Using RSRG with TM method we unravel a 
family of extended (resonant) eigenstates even through 
the underlying lattice lacks translational order. It is not 
unusual that, by controlling the length of a stub appro- 
priately one can have an infinite number of wave- vectors 
of the propagating electron which fall in the spectral gaps 
of an infinite CSW. However, a very interesting feature 
of a CSW array is that, the same wave- vectors corre- 
spond to perfectly allowed (in some cases, even periodic) 
wave- functions for a CSW of arbitrarily large, but fi- 
nite length clamped between semi-infinite ordered leads. 
The most remarkable character of such states is that, the 
nearest neighbour hopping integral on the mapped lattice 
[Fig. 1(b)] exhibit a two-cycle fixed point under RSRG op- 
eration, indicating an 'extended' character (in the RSRG 
sense), and yet the transmission coefficient of the 'lead- 
sample-lead' system decays following a power law, as the 
sample grows in size. We term such states as 'atypically 
extended'. This result is in total contrast to recent stud- 
ies on a Fibonacci quantum waveguide network |l9l l20l | . 

In what follows we describe our work. In section II, we 
introduce the model and the method, section III consists 
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FIG. 1: (a) Second generation Cantor array of quantum stub 
clamped in between two semi-infinite leads, (b) The equiva- 
lent one-dimensional chain. 



of a discussion of the overall transmission characteris- 
tics with an emphasis on the resonant eigenstates and 
in section IV we talk about the 'atypical' eigenfunctions 
and the possibility of the trapping an excitation in such 
a quantum waveguide network. In section V we draw 
conclusion. 



II. THE MODEL AND THE METHOD 

A triadic Cantor sequence is a binary arrangement of 
two elements, A and B for example [Tj|. The sequence 
grows in a recursive manner, starting with a seed A, us- 
ing a substitution rule, A ABA and B BBB. 
The first few generations are. Go = A, Gi — ABA, 
G2 = ABABBBABA, and so on. It is to be noted that 
such a sequence differs from conventional binary quasi- 
periodic sequence in the sense that, in the thermo- 
dynamic limit the sub-clusters of B will span the entire 
one-dimensional space with the elements A acting as im- 
purities located at specific sites determined by the se- 
quence. 

We model our system by attaching a series of stubs per- 
pendicularly to a one dimensional wire. Stubs of lengths 
La and Lb are arranged in a Cantor sequence with inter- 
stub spacings dA and ds as shown in Fig. 1(a). We con- 
sider an electron entering the system from one side of the 
nano-wire with wave-vector q. The wave-function in the 
nth segment (chosen to be the segment between the nth 
and the (n + l)th stubs along the X-direction) is given 
by, 

^„(x) = ^„e''(^-^") + S^e-'"'^^-^") (1) 

The wave-function in the nth stub (along the In- 
direction) is written as, 

My) = C„e''?(^-^") + A.e-^'^^^-^"^ (2) 

The above equations are written using a local coordinate 
system in the spirit of Ref. |JL^. It is then simple to 
match the wave-functions and their derivatives at the 
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nodes to arrive at the matrix equation: 



An+i 1 _ r) I 

Bn+l j " I Br, 



(3) 



where, 
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(4) 



Instead of deaHng with the matrices Q we discretize 
Eq. (3) into a set of difference equation much in the spirit 
of Poincare mapping [J. B. Sokoloff, in Ref. The 
resulting set of equations assume a rather simple form, 



{E 



t 



n,n+ 



n.7i— 1 ^ n — 1 



(5) 



where, ^„ is the amplitude of the wave function at the 
node n. This allows us to deal with the problem in the 
same manner as we do for a tight binding electron of 
energy E moving in a one dimensional lattice with on- 
site potential e„ for the nth site and the nearest neigh- 
bour hopping integrals tn,n±i- In our case the quantity 
{E — Cn) assumes three different values depending on the 
nearest neighbour environment. The values are, 



E~e^ 



2 cot qdB + cot qLb 

cot qdA + cot qds + cot qLa 

cot qd A + cot qds + cot qLg (6) 



The nearest neighbour hopping integrals on the mapped 
lattice are given hy Ia — ^/ sin qdA and Ib = 1/ sin qds 
across the lengths dA and ds respectively. Eq. (15) is 
easily cast into the matrix form, viz. 
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2: Variation of T with q when dA = ds ~ a = 1, 
0.5, Lb = 1. 



TRANSMISSION COEFFICIENT AND THE 
RESONANT EIGENSTATES 
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where we define the three 2x2 transfer matrices M„ as, 
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according to the nature of the nth site (i.e. a, (3, 7). With 
the parameters e/3, and tA, ts and the transfer 
matrices Ma, Mp, and we are now in a position to 
analyze the electron states and transmission properties 
of the equivalent one-dimensional lattice [Fig. 1(b)] using 
the RSRG method. 



We now present the results of our calculation in rela- 
tion to the transmission properties of finite CSW of ar- 
bitrary lengths. Let us, without losing any physics, stick 
to a particular model of the waveguide array, consisting 
of equispaced stubs of two different lengths La and L b ■ 
The equivalent one dimensional lattice then consists of 
two different on-site potentials ca and €3, and a constant 
hopping integral t. This geometry is easily achieved by 
setting — — tB, = and dA = dB = a so that, 
tA ^ tB = t. The transfer matrices are now of two types, 
viz. Ala = = Ma and Mp = Mb- The product 
(7) transfer matrix for a CSW at the ^th generation is given 
by Ml = M/_iM|'"'M;_i, with Mi = MaMbMa- 

In order to evaluate the transmission coefficient of a Can- 
tor stub waveguide (CSW) at any generation, we place 
the equivalent one-dimensional sample between two semi- 
infinite leads. The lead, in a tight binding description, is 
characterized by identical on-site potentials eo and con- 
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stant hopping integral tg. The transmission coefficient 
for such a system then can be worked out in terms of the 
traces and anti-traces of the relevant transfer matri- 
ces and for a Ith generation Cantor stub array it is given 
by, 



T, 



4 sin qa 



[zi cos qa ~ yiY + xf sin qa 



(8) 



where, xi = TrMi, yi = Mi(2,l) - Mi(l,2) and zi = 
A//(l, 1)— M/(2, 2) respectively. The trace and anti-traces 
can be evaluated recursively by using the equations, 



xi+i ^ Ufi-i{xb)[xiwi 

yi+i = UF,-i{xB)[yiwi 

Zl+l = UFi-lixB)[ziWi - 

Wi+i = UFi-l{xB)[wf ~ 



- Xb] - UFi-2{xB)[xf - 2] 

f ys] - UFi-2{xB)xiyi 
\- zb] - UFi-2{xB)xiyi 
2] - Ufi-2{xb)[xiwi - Xb] 



(9) 



Here, = 3' ^ and Um is the mth order Cheby- 
shev polynomial of the second kind, xb — Tr{MB), 
wi = Tr{MBMi), yB = Mb{2, 1) - Ms(l,2), and zb = 
Mb(1,1) — Mb(2, 2) respectively. The initial values of 
the traces and anti-traces are 
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In Fig. 2 we present the variation of the transmission co- 
efficient against the wave vector ranging from zero to 
271 /a for the second and fourth generation Cantor stub 
waveguide system. It is seen that there is a wide central 
gap is the spectrum with the finite transmission zones 
clustering around qa = ^ and With increasing size, 
the two sub-bands around ^ and ^ start getting frag- 
mented [Fig. 2(b)], giving rise to multiple closely spaced 
zones of high transmission separated from each other by 
gaps. Keeping all other parameters fixed as in Fig. 2, a 
change in the value of one of the stub lengths, say. La, 
is found to produce notable changes in the transmission 
spectrum as well as in the nature of the eigenfunction. 
This aspect may be used in a possible experimental setup 
in order to control the transport at a specified value of the 
wave- vector. Let us clarify by giving an example. We set 
U2{xb) = 0. The corresponding value of the wave- vector 
satisfies the equation 

cos qa ~ ± — 
3 
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FIG. 3: T vs La with q = cos^^( — |), dA = ds — a = 1, and 
Lb = 1 (a) second generation, (b) third generation (c) fourth 
generation. 



It can easily be checked that the above choice of q (or qa) 
makes Mg = ±1. This implies that an arbitrarily large 
Cantor sequence of stubs of La and Lb becomes indis- 
tinguishable from a periodic array of the triplet LaLbLa 
If q = i cos~^(±i) happens to be in the allowed energy 
spectrum of an infinite periodic chain of the triplet, we 
get an extended eigenfunction, even though the original 
Cantor sequence lacks translational invariance. In Fig. 3 
we plot the variation of T with La for q — cos~^(— |) 
and a = 1, Lb = 1 within half a period. An interesting 
feature in Fig. 3 is that for such a value of the wave- vector 
to correspond to a finite appreciable value of the trans- 
mission coefficient. La has to exceed a threshold value. 
We have not been able to determine the threshold ex- 
actly, but careful numerical search has revealed that, as 
one deals with longer and longer generation segments, 
the threshold remains close to g = 0.2 (in the unit of a). 
For bigger generations change in the threshold value of 
La becomes practically un-noticiblc 

It is instructive to look at the distribution of amplitudes 
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tials and the hopping integrals are found to be, 
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FIG. 4: Lattice-like profile of an extended eigenstate where 
q — cos~^(|). We have taken dA ~ ds = a — 1, La ~ 0.5 
and Lb = 1. 



of the wave function for the above value of q for a semi- 
infinite stub array. It is well known that the clusters of 
impurities in a one dimensional chain are likely to give 
rise to unscattered (resonant) eigenstates. This is true 
for a random dimer model j22j and its extensions, and 
for a certain class of quasi-periodic chains as well ■ In 
the same spirit one can search for g- values corresponding 
to extended wave-function in a CSW network by mak- 
ing higher powers of Mr equal to identity matrix. The 
method is well known 22, 2^, and we do not repeat it to 
save space. In the example discussed above, the ampli- 
tude profile reflects a perfectly extended character, and 
assumes an interesting lattice-like distribution for a 
specific value of La = 0.5. In Fig. 4 we show such a 
lattice-like profile of an extended eigenstate. Other com- 
binations of the wave-vector and the system parameters 
also give rise to similar observations. However whether 
the amplitude distribution resemble the lattice or not de- 
pends on the specific values of the parameters, the wave- 
function may be of an extended nature. 



IV. ATYPICAL "EXTENDED STATES" AND 
ANOMALOUS TRANSPORT 
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where, 

It is well known that a flow of the hopping integrals under 
successive RSRG steps provide information about the na- 
ture of the wave-function. For example, for an extended 
wave-function the hopping integrals never flow to zero 
under any number of RSRG iterations whereas, for a lo- 
calized state at least one of the two hopping integral will 
flow to zero under RSRG. However, for the latter case 
one has to ensure that the chosen energy does not fall 
into a gap in the spectrum of the infinite lattice. 

We stick to the discussion of a model consisting, as be- 
fore, of two different equispaced stubs of lengths L^ and 
Lb- An RSRG treatment of this model needs the general 
recursion relations (10) described before. The present 
model can be retrieved from the general model described 
before by setting cIa = ds = a and ea = 7^ e/j. We 
now have tA=tB = — . Let us now choose E — en = 

^ ^ sm qa " 

which results in the equation 



cot qLb + 2 cot qa — 



(11) 



This condition is easily be achieved, for a given value 
of the lattice spacing a, by appropriately choosing the 
stub length Lb- For example, with a = 1 and qa = 
{2n + 1)tt/2, n = 0, 1, 2, Lb can assume values equal 
to 2n+i (i'^ unit of a) for which the above equation is sat- 
isfied. The length of the other stub La, of course remains 
free to be chosen. Even though we start up with the spe- 
cial model, on decimation, we land up with a mixed one 
and one has to use the set of recursion relation (10) for 
drawing conclusion. Iterating the set of Eq.(lO) we ob- 
serve that and remain fixed at their initial values 
while iyi(= ^b) exhibit a two cycle fixed point behaviour, 
viz, 

Ia — -tA — Ia 



We begin this section by noting that an infinite Cantor 
stub array can be renormalized into a scaled version of 
itself simply by using the growth rule in the opposite di- 
rection. In terms of the equivalent one-dimensional chain 
this implies the decimation of a subset of sites. The re- 
sulting RSRG recursion relations for the on-site poten- 



tB = -t-B = tB 

Most interestingly increases on successive renormal- 
ization following the prescription 

e^(n) = 2ep{n - 1) - 3e„, 
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FIG. 5: Amplitude ^'jv growth pattern for (a) La = 0.25, (b) 
La = 0.5 and (c) La = .75 for g = ^ (in unit of inter-lattice 
spacing). Here, La = (Ia = (Ib = 1- 



n indicating the stage of renormalization. This impUes 
that an electron propagating in the Cantor stub network 
with above specification will experience, in a renormaliza- 
tion group sense, progressively higher potential barriers 
resulting in a decay in the value of the transmission coef- 
ficient as the system grows in size. A pertinent question 
in this regard is, whether q = (2n + l)7r/2 belongs to the 
spectrum of an infinite or semi-infinite chain. The an- 
swer is 'no'. With respect to the amplitude on the first 
site, the amplitudes at distant sites keep on increasing 
for larger and larger system and consistent solutions of 
the difference equation (5) is not possible to obtain in 
this case. However, it is interesting to note, and simple 
to check that, the same set of all the parameters yields a 
perfectly allowed, even a periodic distribution of ijji when 
any arbitrary generation of a Cantor segment is clamped 
between two semi-infinite perfect leads, where the leads 
are described by identical scatterers placed periodically. 
The growth of the potential under RSRG implies that 
an electron released in the middle of a large cluster of B 
nodes will experience high potential walls on either side 




FIG. 6: Variation of the maximum value of the amplitude in 

a Cantor waveguide as a function of the stub length La- The 
solid and the dashed curves correspond to (a.) q = ^ and (b) 
q = Other parameters arc Lb = 1, dA = da = 1- We 
show I ^'242 i up to 150 for clarity. 

and will be ultimately trapped in this well if the system 
size is large enough. This fact co-exists with the fixed 
point behaviour of the nearest neighbour hopping inte- 
grals which usually characterize an extended eigenstate. 
We thus come across a set of, what we may call, 'atypical' 
eigenstates. 

Both the growth of the amplitude and the decaying trans- 
mission coefficient can easily be estimated if we note that 
for q = {2n + l)n/2, and cIa = ds — Lb = 1, we get 
= —I, so that any arbitrarily large generation of a 
cantor segment of A's and B's become identical, from 
the standpoint of the propagating electron, to a periodic 
array of the unit cells 'AB\ the transfer matrix across a 
finite Ith generation array now reads 

^'"V ) 

where, Nab{1) is the number of AB-h\ocks left in the Ith 
generation after the central clusters of B yielding iden- 
tity matrices. This leads to the power law decay in the 
transmission coefficient as a function of Nab{1), viz. 

Before ending this section, it is worth mentioning that, 

by timing the length of the other stub La one can reduce 
the growth of the wave-functions for g-values which fall 
in the spectral gap of a semi-infinite Cantor sequence of 
stubs. The amplitude of such a wave- function usually 
exhibits a slowly increasing pattern interspaced with flat 
constant regions corresponding to the periodic clusters of 
_B's in between the A's. The maximum amplitude for q = 
^, say, at any generation I occurs at the (3' — l)th site as 
can ciasily be checked. In Fig. 5 we exhibit the flattening 
of the amplitude growth pattern for three different values 
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FIG. 7: Variation of transmission coefficient Eis a function of 
La for q = ^. Wo have taken Lb = 1, = ds = 1- Solid, 
dashed, broken and dotted curve correspond to first, second, 
third and forth generation CSW respectively. 

of La- By changing La continuously from to 1, wc 
get a smooth decay of the largest value of the amplitude 
occurring in any generation until with La = Lb = 1, we 
get a completely periodic distribution of . At this stage 
the lattice becomes indistinguishable from a periodic one. 
Fig. 6 shows how the largest value of the amplitude in a 
243-stub array drops from large values to ±1. 



The same aspect is reflected in the variation of T as a 
function of the stub-length {La, for example). Fig. 7 
shows that perfect resonance occurs at La = 1 only which 
represents a periodic array of stubs. For all other values 
of La the transmittance decays as the system grows in 
size (i.e. as Nab increases). The 'switch-over' to com- 
plete transparency occurs sharply at La = 1 only when 
the system is infinitely large i.e. an 'off' to 'on' state can 
be obtained as La is varied across the value unity. 



V. CONCLUSION 

In conclusion, we have used the renormalization group 
idea and a transfer matrix method, to examine the elec- 
tronic states and transport in finite and infinite stubbed 
waveguide network constructed following a triadic Can- 
tor geometry. Apart from usual extended eigenfunc- 
tions which correspond to complete transmission, such a 
waveguide network when clamped between two perfectly 
ordered semi-infinite leads, is found to support 'atypical' 
states which may display a periodic pattern of amplitude 
and yet the end-to-end transmission decays as the sys- 
tem grows in size. We provide our explanation from a 
renormalization group point of view. 
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